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Abstract 

The formalism of exact ID quantization is reviewed in detail and applied to the 
spectral study of three concrete Schrodinger Hamiltonians [— d 2 /dg 2 + l/(g)] ± on the 
half-line {q > 0}, with a Dirichlet ( — ) or Neumann (+) condition at q = 0. Em- 
phasis is put on the analytical investigation of the spectral determinants and spectral 
zeta functions with respect to singular perturbation parameters. We first discuss the 
homogeneous potential V(q) = q N as N — > +00 vs its (solvable) N = 00 limit (an 
infinite square well): useful distinctions are established between regular and singu- 
lar behaviours of spectral quantities; various identities among the square-well spectral 
functions are unraveled as limits of finite- A properties. The second model is the quartic 
anharmonic oscillator: its zero-energy spectral determinants det(— d 2 /dq 2 + q 4 + vq 2 ) ± 
are explicitly analyzed in detail, revealing many special values, algebraic identities 
between Taylor coefficients, and functional equations of a quartic type coupled to 
asymptotic v — > +00 properties of Airy type. The third study addresses the potentials 
V(q) = q N + vq N ^ 2 ~ 1 of even degree: their zero-energy spectral determinants prove 
computable in closed form, and the generalized eigenvalue problems with v as spec- 
tral variable admit exact quantization formulae which are perfect extensions of the 
harmonic oscillator case (corresponding to TV = 2); these results probably reflect the 
presence of supersymmetric 2 potentials in the family above. 

1 Corrections and updates (July 2003) indicated by footnotes. 

2 In the printed version, we mistakenly used "quasi-exactly solvable" for "supersymmetric" 
potentials throughout (the two notions happen to agree for N = 2 and 6, but not beyond). 
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Exact quantization, or exact WKB analysis, supplies new tools for the an- 
alytical study of the ID Schrodingcr equation, now including arbitrary polyno- 
mial potentials. Here we initiate applications of an exact method to miscella- 
neous concrete problems and models of analytical interest, emphasizing exact 
and asymptotic relations for the spectral determinants and related spectral zeta 
functions. 

We have chosen three rather different quantum potentials to illustrate a 
variety of situations. These have some basic common features (besides their 
required one-dimensional and polynomial nature): they are rather simple, with 
a minimal number of parameters, to remain concretely manageable; one crucial 
parameter (discrete or continuous) governs the transition to a singular limit, cre- 
ating an interesting dynamical and analytical situation; some uniform principles 
for tackling those problems can be issued at the most general level. 

In contrast to earlier studies concerned with individual eigenvalue or eigen- 
function behaviour, we seek the limiting properties of spectral functions, which 
are symmetric functions of all the eigenvalues at once. We take a semi-rigorous 
approach, in which we argue a global operational scheme without claiming ab- 
solute completeness in every detail. 

The paper is organized as follows. 

The introductory Sec.l gives a detailed survey of the exact tools to be used 
here, while an Appendix collects all the results concerning homogeneous po- 
tentials. This is done partly for convenience, given the lack of comprehensive 
reviews for results that are very scattered in time and place of publication, 
and also to clarify some parts of our most recent developments where localized 
inconsistencies went undetected. 

Then, Sec. 2 (Exercise I) considers the family of homogeneous potentials 
— \l\ N as t ne degree N tends to +00. This is a most singular problem, 
for which many explicit results are however available beforehand, and the lim- 
iting problem (an infinite square well) is exactly solved by elementary means. 
We therefore mainly propose and test some general principles of investigation, 
rather than claim truly new results. In particular, we suggest criteria for sort- 
ing out regular vs singular types of limiting behaviour in spectral zeta functions 
and determinants. Still, we identify several possibly unnoticed properties and 
formulae in the N = 00 limit which arise as regular limits of nontrivial finite- N 
properties. 

Sec. 3 (Exercise II) deals with the quartic anharmonic oscillator family V(q) = 
q 4 + vq 2 which is the most common model for singular perturbation theory (the 
free harmonic oscillator emerges in the v — > +00 limit). We single out a pair 
of one-parameter spectral functions for their remarkably numerous and sim- 
ple explicit properties: the zero-energy determinants Qi ± (i>) = f det(— d 2 /dq 2 + 
q 4 + vq 2 ^ (+ corresponds to the even-state sector, — to the odd-state sec- 
tor) . We present a simple WKB technique allowing to express asymptotic rela- 
tions between ^-dependent determinants such as det(— d 2 /dq 2 + q 4 + vq 2 )* 1 and 
dct(— d 2 /dq 2 + vq 2 ^ when v — > +00. Then, practically all the analytical re- 
sults available for the homogeneous quartic case (V(q) = q 4 ) have counterparts 
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for the functions Qi (and their associated spectral zeta functions), while their 
asymptotic properties remind of the Airy functions (see fig.l). Several special 
values are computable and listed in Table 1 against the analogous results for 
the Airy functions and the quartic determinants det(— d 2 /dq 2 + q 4 + X) ■ 

The final Sec. 4 (Exercise III) gives a complete treatment of the similar deter- 
minants for a different class of binomial potentials, V(q) = q N + vq N / 2 ~ x (for TV 
even) . Here the formalism yields fully closed forms for the zero-energy spectral 
determinants: ea. 1)121)1) for det[— d 2 /dq 2 + V^dgl)]* in terms of Gamma func- 
tions (plus non-trivial exponential prefactors) and, when N is multiple of 4, the 
still simpler eq. l|123[) for the zero-energy determinant of the potential V(q) itself 
on the whole real line; exact quantization formulae follow for the correspond- 
ing generalized spectra (in the v variable: eas. (|107|l and (|124fl respectively). 
A broad generalization of the familiar harmonic-oscillator exact results is thus 
obtained; for N = 2 [mod 4], this seems to describe a zero-energy cross-section 
of the formalism for supersymmetric 2 (including some quasi-exactly solvable) 
models. 

Although all three problems rely on the same background formalism, they 
can be approached fairly independently from one another. Accordingly, there 
are no global conclusions but each Section carries its own concluding remarks. 

1 Introduction 

1.1 General results on spectral functions [TJ, [2] 

Here, an admissible spectrum is a purely discrete countable set {Ek}k=o, 1,2, ••• 
with Ek > 0, I +00, such that its partition function 

00 

6(t) d = e" tBfe Re(t) > (1) 

k=0 

can be asymptotically expanded in increasing (real) powers {t p } for t \ 0: 

0{t) ~ c P tP (* J- °)' with ^ =' ~ min > ( thc "g rowth order"). (2) 

p 

(In summations etc., Latin indices will systematically mean integers, and Greek 
indices mean more general real indices, namely real-valued functions over the 
natural integers, strictly f +00 or J, —00.) 

1.1.1 Spectral zeta functions [3 |H gj 

Eq.J2J) implies that the "Hurwitz", resp. "plain" spectral zeta functions 
Z(s,\) d ^J2(Ek + \)- s (|arg(A+£ )| <%-6), resp. Z{s) A ^ Z(s, 0) (3) 

k 
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converge for Re(s) > p, that Z(s) has a meromorphic continuation to all com- 
plex s, with 5 

polar set: {— p}, residue formula: lim (s + p)Z(s) = c„/ T(—p), (4) 

s—-p 

"trace identities": for m e N, Z(—m) = (-l) m m! c m (c m = f for m £ {p}) 

(5) 

and similarly for Z(s, A), just by substituting e xt 0(t) for for general A, 
the leading trace identity (specially useful for ea. l|15|) below) is then 

Z{Q,X)= V c _„i^ = Z(0)+ V Res s= „Z(.s) (6) 
1.1.2 Spectral determinants [T] 

Since i/ (s, A) is regular at s = 0, a spectral determinant can be defined by 
zeta-regularization, as 

D{\) = det(# + A) =' cxp[-a s Z( s , A)] s=0 ; (7) 

it is an entire function of order p with {— Ek} as its set of zeros. Moreover, 
amidst all such functions, D(X) can be precisely picked out in at least two ways. 

• On one side, eq.J2J implies a canonical semiclassical behaviour for -D(A), 

-logD(A) ~ ^ Cp r p (A)forA^+oo, r„(A) d = d s [ ^±^ A~^ 

1 - 6 -' iplAj ~ l(-(-A) m /m!)(logA-Er=i 1 A) if-p = meN W 

are the only terms allowed; no other type, including additive constants (oc A ), 
can enter this expansion. 

• Independently, D(X) is also fully specified by expansions around A = 0: 
firstly, in reference to the Fredholm determinant A(A) (built as a Weierstrass 
infinite product): 

D(X) = exp[-Z'(0)- V ^(-A)"lA(A), (9) 

l<n</x 

oc 



where A(A) d ^ fl (l + ^) cxp [ £ (for all A), (10) 

fc=0 k l<n<fj, k 



and Z{n) — Z(n) if Z(s) is regular at n, as when n > p (11) 

Z(l) = lim (Z(s) - — (= Z{1) or its finite part) (12) 
s— *i V s — 1 / 

(the general formula for Z(n) on a pole is more contrived and required only 
for p > 2 (P], eq.(4.12)), whereas p < 3/2 in this work). Lastly, by way 
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of consequence, the determinants are also characterized by these Taylor series 
(converging for |A| < E ), 

-log£(A) = Z'(0) + f^(-A)", -logA(A) = (13) 
* — ' n * — ' n 

n—l n>fi 

The simplest case is fi < 1: then 

A(A) = n( 1 + ^)= cx p[-E^r(- A ) n ]' DW 55 e ~ z ' m A(A) - (14) 

fc=0 fe 71=1 n 



Two other properties are worth mentioning: 

- if a spectrum {Ek} is dilated to {aE^} (a > 0), the spectral functions get 
rescaled to 

Z(s,X\a) =a- s Z(s,\/a) => D(X\a) = a z( -°' x/a) D(X/a) (15) 

(a behaviour hence mainly governed by the leading trace identity, eq.@); 

- finally, all those results extend to analogous complex spectra jH]- 

1.2 ID Schrodinger operators with polynomial potentials 

We subsequently specialize to ID Schrodinger equations involving a polynomial 
potential V(q) (adjusted to V(0) = 0), 00 

(-d 2 /d(? 2 + [V(q) + \})4> = 0, V{q) = +q N + [lower order terms]. (16) 

We call H + (resp. H~) the Schrodinger operator on the half-line {q > 0} 
with the Neumann (resp. Dirichlet) boundary condition at q — 0, and H the 
Schrodinger operator on the whole line with the potential V(M)j whose spec- 
trum we denote by {Ek}- Then {Ek}k even (resp. {Ek}k odd) is the spectrum of 
H + (resp. H~) 1 each one is an admissible spectrum in the previous sense, with 

growth order fx = h + (m < 1 generically, i.e. for N > 2), (17) 

exponents {p} = {-fx + j/N} j=0 ,i, 2 ,-; (18) 

their spectral functions Z + , D + (resp. Z~ , D~) arc our basic concern, and 
have properties as above (with exceptions in the singular case N — 2). However, 
a few results take a neater or a more regular form upon recombined functions 
instead, 

Z d =Z+ + Z~, D d = D+D~ (spectral functions of H), (19) 
z p def z + _ z -^ d p 4ef D +j D - (" skew » spectral functions). (20) 
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1.2.1 Classical "spectral" functions [9j 

The quantum spectral functions of the problem ijPol) admit natural classical 
counterparts with parallel properties. The Weyl-Wigner correspondence, for 
instance, associates the following classical partition function to the quantum 
one of eq.Q, 

Mt ) = /" dgdg e ^ +vm]t _ 1 [ + °° ^ v(q)t d ^ (21) 



2tt y/ni Jo 

with an expansion coinciding with eq.(|2J) as long as p < 0. The same Mellin 
transforms as from the quantum 6{t) to Z(s, A) then yield 

Z cl (s,X) d = f -L [ + °°Ut)e- Xt t s - 1 dt= r ^' 1 /- 2) /o( g ,A)(22) 

r(sj 7 us)^ 

/•+CO 

where I q (s,X) = / (F(g') + A)^ +1/2 dg' (Re(s) > M ); (23) 

and £>d(A) d = exp[-d s Z cl (s, A)] s=0) (24) 

with properties induced by eq.@ similar to the quantum case (but D c i(A) is 
not an entire function: branch cuts replace the discrete zeros of D(X)). 

The meromorphic continuation of I q (s, A) is thus important at s = 0. If (for 
A > — inf V initially) we compute the expansion 

(V(q) + X)- s+1/2 ~ Yl Ms)^ for«-+oo (a = f , ^ - 1, • • •) 



(the (3 a of course also depend on A and on the potential), then 



(25) 



(T 

hence it satisfies 

lim sl q (s, A) = 0- 1 (O)/N = -Z(0,X)/2; (27) 

(3-i{s) is actually independent of A except for N — 2; the latter value for the 
residue directly comes from ea. (l22H . making another explicit statement of the 
trace identity |JBJ. 
We denote here 

X q {X) d = lim jl 8 (s, A)- ^=^ | (the finite part of 7 g (s, A) at s = 0). (28) 

Now, some of our earlier statements revolving around this quantity [9*lll(J| require 
corrections in the most general setting /3_i(s) ^ 0. Finite parts are also to be 
extracted: first on the expansion l|26|) . giving 

T,(A)~- £ /3(7 (0)i^l-/3_ 1 (0)log 9 +^^ (g-+oo), (29) 
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then on the definition I|22I24[) . giving 

~ \ogD cl (X) d ^ -\d s Z cl {s, X) s=0 = I (X) + 2(1 - log 2)^51 . (30) 

The finite partX (A) of / + °° (V(q) + X)- s+1 / 2 dq at s = is a candidate to define 
a "symbolic" value for the divergent integral / +O ° H\(q)dq, where 

H\(q) = (V(q) + A) 1 / 2 , the classical (forbidden-region) momentum. (31) 

However, such an assignment being conventional, we much prefer the "renor- 
malization" given by ea . (|30|) . which only adds an explicit constant to the finite 
part. So, (more generally) we pose the suggestive notation 

" +0 ° - - 1 rq 0_i(O) 



J Tl x (q')dq' ^ - log D cl (A) - ^ Il x (q')dq' = J q (X) + 2(1 - log 



2)- 



N 
(32) 

A big advantage of eq. I|32l) will be that its A — * +oo expansion has the canonical 
form © (essentially because \ogD c \ behaves similarly to its quantum counter- 
part log-D in this respect, and J Q 9 U\(q')dq' = 0(^f~X) + o(l) is also manifestly 
canonical). 

Thanks to ea. l|29(l . the prescription (|32|l is also directly characterized by its 
large-g asymptotic behaviour, as 

-f oo 

IW)<V = -5 A (q)-/3_ 1 (0)log 9 + C + O (l) (q -> +oo) (33) 

■ ■ 1 < I Jti.s) 

s=0 



where S x (q) V /^(O)-^-, C ^ ^-(-2 log 2/3^(0) + d s 

(T + l TV V 

{<r>-i} 



1 - 2s 

Classical analogs will also arise for D + and D~ separately fea. l|46|l below). 



1.2.2 Special features of ID Schrodinger determinants 9i ll()| 

The quantum determinants D ± for eq. l|16|l can be specified in two additional 
ways. 

• Either spectrum of H + or H~ obeys a semiclassical (high-energy Bohr- 
Sommerfeld) quantization condition of the form 

££,V~* + l/2 for c o v ™fc^+oo, (34) 
p 

which implies Euler-Maclaurin continuation formulae for the zeta functions; 
e.g., down to Re(s) > 0, 



Z ± (s) = lim ( V E- s + -E K S - - V f—^E K s - p \ for ov ™ k, K . 
v ' K ^, +00 \Z-^ k 2 K 2 ^ s + p K J odd ' 

fc<A' {p<0} P 

(35) 
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The resulting residues are compatible with cq.Q provided 

(the latter specifies the s = trace identities); here, moreover, 

bt p = bZ p for all p < =>> Z(0)=6jjS Z p (0) = 1/2; (37) 

then the same trace identity follows for all Z(s, A) and Z p (s, A) by eq.lJBJ, except 
for N — 2 (eq.JZU below). 

By the same token, D ± (X) become directly specifiable as functionals of the 
spectrum: in the simplest case (J, < 1, i.e. N > 2, 

logD ± (\) = lim { V log(£ ft + A) + \ \og{E K + A) (38) 

k<K 



2 

{p<0} 



even 
Ddd 



(the expansion coefficients b_ p being themselves functions of the spectrum {Ek})- 

• D ± (X) are also related to an exact solution ijj\{q) of ea. l(Trj|l defined by a 
particular WKB normalization, 



/ / \ tt / N-l/2 „ f + °° n *0?')dg' 



for n A (g) -> +oo; (39) 



here n>(g) is the classical momentum function l|31|l . and the divergent integral 
f °° U\(q')dq' is specifically "renormalized" through eq. ; however, eq.lO 
also makes ij)\(q) directly specified as the (unique) solution of the Schrodinger 
equation H16J) that decays for q — > +oo with the precise behaviour 

Mq)~e C q~ m ~P- li0) e~ Sx{q H=e C q NZ ~ {0 - X) e- Sx{q) ), q^+oo (40) 
(the latter form comes from eas. l|27l37|l ). Remark: in parallel to eas. (|39l40|l . 

-i,' x {q) ~ n x {q)^J^ nM)Aq ' for U x (q) -> +oo (41) 
( ~ ^ ? +nr/i-fl-i(0) e -5x(g) =e CgiVz+(o,A) e -s A ( g ) for(? ^ +0O ). 

We refer to ip\(q) as the "canonical recessive" solution. (It is proportional 
to the "subdominant" solution of 7 , ch.2, but only equal to it if /3_i(s) = 0. 
The current normalization also differs from and from |!3 EH — where it 
suffers localized inconsistencies. The discrepancies, which for N > 2 involve 
A-independent factors only, ultimately cancel out in all results involving only 
spectral determinants.) 
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The extension o{ip\(q) to the whole real line, as a solution of eci. (|16[) with the 
potential V^(|g|), is shown by integrations to identically satisfy Apps.A,D) 

D(X) = CW X: C p (A) = C p [-^(0)M(0)], (42) 

W x d = Wronskian{^ A (- <Z ),^ A (g)} = -2^(0)^(0) (43) 

for some constants C, C p ; to identify these, we test a characteristic prop- 
erty of the spectral determinants: the canonical large- A asymptotics (eq.©) of 
their logarithms down to constants included; as regards the right-hand sides in 
ca. l|42l) . the WKB formulae 139I41|) (good for large A) supply these asymptotic 
forms, 

WV2~e 2 /o +00n ^ d « = i} cl (A), -V>A(0)M(0)~n A (0) (A^+oo), (44) 

both of which have canonical (A — > +oo) logarithms (D C \(X) by analogy with 
D(X), and n A (0) by inspection); hence necessarily, 

D(A) = W x /2 = -^(0)^(0), D P (X) = -^(0)/Va(0). (45) 

Incidentally, the second ea.l|44|) then naturally specifies a "classical skew deter- 
minant", -D P i(A) = f n A (0), from which classical analogs of D ± (X) also follow 
as 

D±(A) d ^ Il^Of 1 / 2 eX + °° n ^> d<? . (46) 

The main conclusion however concerns the quantum spectral determinants 
D themselves: upon a straightforward simplification of eq. l|45|l . they get ex- 
pressed in terms of ip\ (q) by the fundamental identities 

D-(A)=Va(0), D+(X) = -i// x (0), (47) 

(also valid for a rescaled potential, i.e., V(q) — vq N + ••■)• 

1.2.3 The main functional relation 9l ll(J| 

Jointly with the original problem its set of "conjugate" equations is defined 
by means of complex rotations, as ( 7 , ch.2) 

VW(q) = e-^V(e-^ 2 q), X [e] = f e"** A, <p = ' (48) 

where £ = 0, 1, - ■ ■ ,L — 1 [mod L] labels the distinct conjugates, of total number 

N 

L = N + 2 in general, L = — + 1 for even polynomial potentials. (49) 

The main result to be used throughout is the WronsJcian identity, which 
states a bilinear functional relation between the spectral determinants D ± (X) 
and those of the first conjugate equation, namely D^ ± {e^ lip A): 

e +i V /4 D [l] + (e - i(P A )£,[0]-( A )_ e -i V /4 £) [0] + (A)D [l]- (e -i V A) _ 2ie i^_ 1 (0)/2_ 

(50) 
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It entails an exact quantization formula for the eigenvalues 



1 argDl^C- e~» E) E=Ek - f ^(0) = fc + i ± for fc = JJfc. 

(51) 

i.e., this condition determines the spectrum {Ek} exactly in terms of the spec- 
trum {-E^ 1 '} of the first conjugate potential (of which the left-hand-side 
£)[l]±(_ e ~ l <p E) are functionals). (Eq.JSlJ together with all its conjugates ap- 
pear to form a determined system for the resolution of all the spectra {-E^ } a ^ 
once.) 

Here we will often invoke the specific results relating to the homogeneous 
potentials 

H N = - — + q N , q e [0, +oo) N > 1 integer. (52) 

The corresponding formulae are collated in greater detail separately in Ap- 
pendix, and also in Table 1 below for N = 1, 4. 

1.3 Statement of the problem 

The previous exact analysis of Sec. 1.2 is controlled in an essential manner by 
the degree N of the potential, quite sensibly since q N defines the most singular 
interaction term. Still, there are interesting transitional situations where this 
parameter N can diverge or behave discontinuously while the quantum prob- 
lem itself has a well-defined limit (example: the term q N is multiplied by a 
coupling constant g — * 0). Earlier studies of such problems show those limits 
to be very singular, and non-uniform over the energy range. An unexplored 
challenge is then to control the limiting behaviour of spectral functions, which 
are symmetric functions involving all eigenvalues simultaneously. The earlier 
analyses, which essentially work at fixed N, need further development to handle 
such non-uniform regimes. 

Here we will begin to gather some insight about this issue by examining 
three model problems at varying depths. 

2 Exercise I: infinite square well as N — > +00 
limit 

The infinite square- well potential over the interval [— 1,+1] can be realized as 
the N — > +00 limit of the homogeneous potential V(q) = \q\ . The solutions 
of the Schrodinger equation with this potential indeed approach those of the 
infinite square well, but the limiting behaviours are interestingly singular, and 
non- uniform with respect to the quantum number k [T2j . 
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2.1 The N = +00 problem 

The limiting Schrodinger operator is given by V(q) = in [— 1, +1] with 
Dirichlet boundary conditions at q — ±1. At infinite N, the spectrum of 1^1^ 
becomes explicit again, {E k — (k + l) 2 7r 2 /4}, of growth order fj,^ = ~. As with 
finite N, the operator splits into (over the even eigenfunctions, labeled 
by even k) and H^, (over the odd eigenfunctions, labeled by odd k). 

Here the immediately explicit spectral functions are the Fredholm determi- 
nants A* of ea. (|14fl (they reduce to standard Weierstrass products), 

. _i_/.s / — r . /,n sin v— A . ... sin2V— A 

A+(A)=cos^, A^(A) = — =^> A oc (A)= (53) 

and the plain spectral zeta functions, 

Z+(s) = ^^a^), Z oc (,s) = 4iC(2 S ) z oc .( S )=(-) 2s a2 S ). 

ir zs ir zs \7r/ 

(54) 

The latter formulae imply the explicit computability of Z^(s) at all integers 
s e Z, e.g., 

Z+(0)=Q, Z"(Q) = -l/2 => Z ao (0) = -1/2, (55) 

and ultimately that of the spectral determinants themselves (even though these 
will not serve here) through D^(X) — exp[— (Z^,)'(0)]A^(A), using 

exp[-(Z±) / (0)]= J D±(0)=2 => £>«,(()) =4. (56) 

2.2 The transitional behaviour problem 

Regarding the behaviour of the spectral functions, a first task is to seek condi- 
tions ensuring the regular behaviour of a quantity (meaning that it has a finite 
limit, and this is the correct value for the limiting problem). For the other (sin- 
gular or pathological) quantities then comes the additional task of describing 
their precise behaviours. 

A fundamental quantity in these problems is the growth order [i. Here, the 
limiting (square- well) value /Uqo = \ agrees with the limit of [j,n (= \ + jf by 
eq.^TJ); hence fj, behaves regularly (as opposed to the later examples). 

As a crude dividing line between (generic) singular and regular behaviours, 
we expect that essentially those quantities which converge both for finite N and 
in the limiting problem should be regular, in particular: 

Z±(s,X) ^Z±(s,X) iff Re(s) > (57) 

Singular behaviour should then set in at Re(s) = and, plausibly, become 
worse as Re(s) decreases further. This purely qualitative argument cannot how- 
ever predict the precise behaviour of any singular quantity. 

Fortunately, quantitative statements are made easier here by a set of explicit 
results for the finite-A problems (cf. Appendix A.l) and their counterparts 
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for N = +00 (Sec. 2.1 above). Those data indeed behave consistently with 
the prediction (|57|l taken with fi^, — \- As examples of regular (s > /ioo) 
behaviours, Z^(l) — > Z^(l) — j/g , and likewise for the higher-order sum 
rules involving s = 2, 3, • • • (Sec. 2. 5 below). As opposite examples (involving 
s < fioo): Z N {0) = for all N whereas Z^Q) = -1/2 by eas J128l55|) : and 
worse, ea. 1(136(1 implies 

exp[-(Z±)'( S = 0)] = D%(\ = 0) ~ (N/tt) 1 / 2 00, (58) 

even though 1)^,(0) have finite values, perfectly defined by ea. (|56(l . in the N = 
+00 problem ! (All that shows how carefully such transitional problems must 
be handled. — The same formulae show the skew spectral functions 1)20(1 to 
behave slightly better: e.g., Z&(0) = 1/2 = Z£(0), D P N (0) - 1 = l£(0).) 
We then basically expect 

A±(A)^A±(A) but D^(X) diverge, (59) 

this divergence being confined here to the factor exp[— Z^(0)] = £>jv(0) alone, 
because fi^ < 1 and the Fredholm determinants can be expressed using s = 
1,2, •• • only as in cq. (|14l> . (The Fredholm determinants A(A) should behave 
regularly in general since they are designed by retaining only regular values 
Z(s) (at integers s) in their Taylor series (cf. eg . lfT3"ll ) : this has to be qualified 
only if fx reaches (or jumps across) an integer in the limit, as in the example of 
Sec.3.) 

2.3 The main functional relation 

We now study the N — * +00 limit of the functional relation (|50|l specialized to 
homogeneous potentials \q\ , as stated in ea. l(13Q|l . By the preceding arguments, 
this functional relation should be well-behaved as N — > +00 only once it has 
been transcribed for Fredholm determinants (using eas. fl4ll?t2f> ): 

e+^/4 A +( e -iv A ) A -(A)-e-^/ 4 A^(A)A^(e-^A) = 2isin <^/4. (60) 

In the latter formula, holding for all ip = 4tt/(N + 2), the expansions of both 
sides in powers of ip — > should then be identified order by order. Ea. l(6()|l 
having the form of a "quantum Wronskian" identity for finite <p ^1 El> it 
is not surprising that the identification to the leading order 0(ip) discloses a 
"classical" Wronskian structure: 

A + (»ajA-)-(A;ixA + )A-.i(l-A+A-). (») 

However, we know neither how to interpret the right-hand side, nor how to 
solve this functional relation directly (and identification at the next order in 
tp oc 1/N in ea. l|60|l does not appear to yield any new information either): this 
constitutes an interesting open problem, since the finite-TV equation admits of 
constructive solutions by an exact quantization method using ea. ((133(1 : this 
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method itself however seems totally singular in the N — ► +00 limit. At the 
same time, the Fredholm determinants of the infinite square well are known, 
given by ea. (|53[l : they explicitly verify ea. (|61[l . and this provides a positive test 
of regular N — > +00 behaviour for the main functional relation in the form J5DJ). 



2.4 The N = 00 coboundary and cocycle identities 

For the homogeneous finite-A problem, a closed functional equation for the 
complete determinant D(A) is supplied in Appendix as "the cocycle identity" 
11341) . a sum of L = O(N) terms. Its naive A — ► +00 limit will be an integral 
relation for -D(A), further reducible by the residue calculus. It is however simpler 
to work out the A — > +00 limit directly upon the (logarithm of the) underly- 
ing "coboundary identity" l|131fl equivalent to ea. (j6U|) : this limit is manifestly 
equivalent to eg . 1|> and has the (additive) coboundary form, 

- 24 logA P (A) m 

A p being meromorphic and A entire, the main solvability condition for ea. (|62[l 
is that the residues at the poles of 1/ A must match the explicit residues of the 
left-hand side, resulting in a curious constraint upon A alone at its zeros, 

A-^A(A)1 = (-l) k for all k e N, (63) 

dA lx=-E k 

which stands as N = 00 counterpart for the cocycle identity (|134fl . It is veri- 
fied by Aqo but, as with ea. (|61|l before, we have no idea about other possible 
solutions. 



2.5 The N = 00 sum rules 

The Taylor series of both sides of ea. (|62[) can be expressed with the help of 
ea.lfT3)l. giving 

n— 1 m— 1 r— 1 m— 1 

(64) 

As in the finite-A^ case (ea. p41|l ). this acts as a generating identity: the iden- 
tification of each power A™ in ea . l|64|l yields a sum rule of order n, which here 
expresses the combination 2Z p (n) — Z(ri)/n in terms of the lower Z(m), as 

2Z p (l) - Z(l) = 0, 2Z P (2) - Z{2)/2 = Z{l) 2 /2, (etc.). (65) 

The regular behaviour of eq. I|60(l implies that these sum rules too must be the 
limits of their finite- A counterparts (|142|l and also be verified by Z^(s). (Due 
to the special form (|54|l of Z^(s), these rules amount to equating each Bernoulli 
number Bm to a certain polynomial in its predecessors.) 
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In conclusion, this case provides a testing ground for ideas and methods 
applicable to transitional regimes, however it has not yet yielded any new results 
about the underlying spectral problems themselves. Still, we have identified 
several novel structures in the N = oo problem as imprints of nontrivial finite- 
TV features in the N — > oo limit. There remains to effectively handle the finite- A~ 
problem as a regular deformation of this N — oo case, but this would probably 
require answering the various questions we left open. 

3 Exercise II: anharmonic perturbation theory 
as N = 2 limit 

We now study the approach towards the other singular limit of the formalism, 
N = 2. It cannot be realized through homogeneous polynomials, but a transition 
from N = 4 to N = 2 precisely underlies the well known perturbation theory 
for the anharmonic potentials |15HlbllT7lll8j 

U g (q) = q 2 + gq 4 (g — > + ), or equivalently V v (q) — q A + vq 2 (v — > +00); 

(66) 

we recall the basic unitary equivalence between the two operators 

H d = -d 2 /dq 2 + V v and y^{-d 2 /dq 2 + U g ) with g = v~ 3/2 . (67) 
We denote the ^-dependent spectral functions of ff ± as 

Z ± (s,\;v), Z ± {s;v) d = Z ± (s,0;v), D ± (X;v). (68) 

The poles of these zeta functions lie at s = 3/4, 1/4, —1/4, • • •; moreover, 
by a straightforward computation of ea. (|25J) . /3_x(s) = for an even quartic 
potential; hence the leading trace identities are u-independent, as 

Z(0, A; v) = 0, Z p (0, A; v) = 1/2 for all finite v. (69) 
3.1 The transition v — » 00 

This transition is now discontinuous: N retains the fixed value 4 for all finite v 
while it has the (more singular) value N = 2 at g = 0. All related parameters 
are then singular, especially now the order: /i = | for all finite v, vs /ico = 1. 
Each eigenvalue of H satisfies 

E k ~ Vw(2/c + 1), v -> +00, (70) 

but not uniformly in k |17j , hence it is another matter to find the behaviour of the 
corresponding spectral determinants D (A; v) themselves, as entire functions of 
A and v. 
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Following ea. l(57|) again, we now expect Z ± (s,X;v) to be well-behaved as 
v — ► +00 iff s > fioo = 1. For instance, the s = trace identity behaves 
singularly: for the limiting potential U g= o(q) = q 2 , ea. H59fl gives 



According to ea. lfHIjl . (Z ± )'(0, A; v) and D ± (X;v) should be singular (as previ- 
ously), but now so should Z ± (l,X;v) and the resolvent trace d\ logZ3 ± (A; v) 
(also involving s = 1); thereafter, higher derivatives (d\) n log D ± (A; v) should 
behave regularly (as they only involve s — n, n + 1, • • •). Understanding the 
v — ► +00 behaviour of D (X;v) then just requires the control of two (pairs of) 
functions of v alone, (Z ± )'(0;w) — or equivalently l? ± (0;z;) — and Z ± (l;v). 

From now on we will exclusively deal with D ± (0; v), a problem which entirely 
resides in the {A = 0} plane. It is technically simpler because the anomaly 
l|71|) present in (and only in) the limiting problem U (q) — q 2 vanishes at A = 
0. Moreover, these restricted spectral functions £) ± (Q;t!) will display many 
explicit properties making them intriguingly similar to spectral determinants 
of homogeneous potentials. The analysis of the complete spectral determinant 
D ± (X;v) (including the coefficient Z ± {l]v) as d\ \ogD ± {X\ v)\ =0 ) is also under 
way but will be more involved. 

3.2 "Extraordinairy" spectral functions 

We therefore restrict subsequent attention to the following pair of restricted 
determinants: 



These entire functions of v will display numerous explicit properties. (Determi- 
nants of general binomial potentials, det(— d 2 /dg 2 + <7 W + vq M ), M < N, can 
be handled likewise, cf. Sec. 4.) 

Two immediate results are: a pair of special values (cf. ea. H136(l and Table 



Qi+(0) = D+(0) = Qr(0) = DZ(0) = W%\ (73) 



and a main functional relation, drawn from eq.J5U|) (with /3_x = 0) and from 
the conjugacy formula V™(q) = q 4 + (fv)q 2 (cf. eq.f^SJ)): 



e+i7r /6 Qi+Qt,) Qi-(«) - e"-/ 6 Qi+(w) Qi" Q«) = 2i (j d = e 2i ^ 3 ). (74) 



This functional relation is identical to that of the homogeneous quartic problem 
(|143|) if the even and odd arguments are interchanged, otherwise it also resembles 
the Airy relation l|149|l . Actually, will display hybrid properties between 
those two cases (hence their name). 

3 A copying error made three numerical values for the Airy case (at s = 0) slightly inaccurate 
in the printed version. 



Z(0, A) = -A/2 (but still, Z p (0, A) = 1/2). 



(71) 



Qi± (v) d = D ± (A = 0; v) = dct(-d 2 /d 9 2 + q 4 + vq 2 ) 



(72) 
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Spectra: 


Airy 


zeros 


quartic 


oscillator 


zeros 


of Qi ± 


Z(s): 








zi 


Z+ 


z- 


Z'(0) 


0.0861126 


-0.2299537 


-0.1460318 


-0.5471153 


-0.1685422 


-0.1780313 








2 4 / 3 3 1 / 3 v ^F 




6 1/3 r(i/4) 


x/2r(i/6) 




3V3 r(l/3) 

= D+(0) 
0.9174909 


3 2 /3 r ( 2 /3) 

= DI(0) 
1.2585417 


r(i/6) 

= ^4 + (0) 

1.1572330 


3V3 r (5/6) 

= ^r(o) 

1.7282604 


r(i/6) 
= X>+(0) 
1.1835782 


6 1 /3 r (i/ 4 ) 

= 2?-(0) 
1.1948628 


Z(0) 


1/4 


-1/4 


1/4 


-1/4 


1/8 


-1/8 


Z(l) 







= 2Z 4 "(1) 
1.5266059 


2 4 / 3 7T 3 


3 4/3 r(2/3) 5 


= 2Z+(1) 
-0.3960418 


— T 

— 0.7290111 


3 17/6 r ( 2 /3)5 

0.7633029 


2 10/3 T 2 

-0.1980209 


Z(2) 


1/r 

1.3717212 


r 2 

0.5314572 


cf. ^] App.C, 

0.9147383 


cqs.(C.28,33,34) 

0.0815825 


w 0.3578 


« 0.2377 


Z(3) 


1 


-T 3 + i 

0.1125618 


0.8414950 


0.0190222 


w 0.10338 


w 0.03859 



Table 1: Analytical and numerical zeta- function values for several spectra. 
First two columns: Airy zeros 3 (cf. Appendix A. 2. 2; notation as in ea. l|153|) : 

r d = - Ai'(0)/Ai(0) = 3 1 / 3 r(2/3)/r(l/3)). Middle two columns: quartic os- 
cillator (Hi) levels (cf. Appendix A.2.1). Last two columns: zeros of Qi* (cf. 
Sec. 3. 2; we could only compute Z ± (2) and Z ± (3) numerically from scarce spec- 
tral data, at an unreliable accuracy). Computations used analytical formulae 
when available, but were also cross-checked against direct calculations upon 
numerical spectra. Note: r(l/6) = 2 2 / 3 ^r(l/3)/ r(2/3) = 2 2 / 3 3 1 / 3 Vtt/t. 
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3.2.1 Zeros of Qi* 

The entire functions Qi ± vanish at those values v — —Wk (< 0) for which A = 
is an eigenvalue for the potential q 4 + vq 2 (generalized spectral problem). As 
before, the labeling of u>k in increasing order makes the parity of k match that 
of the eigenfunctions. (We do not expect any complex zeros.) 
The first few zeros of Qi 1 * 1 evaluate as: 

k (Qi+) k (QT) 

-2.2195971 1 -3.2511776 

2 -5.4900693 3 -6.1598396 (75) 

4 -7.9276920 5 -8.4854215 

6 -10.029209 7 -10.525121 

The zeros (— Wk) also obey an exact quantization condition immediately follow- 
ing from ea. ((?4*|) by analogy with its quartic (|144fl and Airy (|150|l counterparts, 
of which it looks as a crossbreed: 

2 11 

- arg Qi* (-jw) w=Wk =k + -±- for fc = (76) 

Anticipating on the next paragraph on asymptotic results, we can at once 
estimate the zeros {wk} for large k by solving the Schrodinger equation with 
potential V v (q) = q 4 + vq 2 and with A = semiclassically, for large v — —w < 0. 
For q^> l,ip(q) has standard WKB forms in both allowed and forbidden regions 
(decaying in the latter); at the same time, for q <C y/w (far inside the allowed 
region), ip(q) must approximately satisfy (— d 2 /d<? 2 — wq 2 )ijj = 0, an equation 
solvable by Bessel functions (cf. ea. H135|) 1: specifically here, 

^±(g)«q 1/2 J Tl/4 (V^g 2 /2) Cdd solutions). (77) 

Matching the two approximations in the intermediate region {1 <C q <^ \/w} 
then yields the semiclassical quantization formula 

If 3 1 

— <£ pdq ~ k H — for k even, k H — for k odd. (78) 

2?r J p 2 +V ( q ) =0 4 4 

(Remark: this Bohr-Sommerfeld quantization rule is appropriate for a general 
symmetric double- well potential with the barrier top precisely kept at the energy 
0; it induces a splitting between even and odd quantized actions exactly half- 
way between their (quasi)degeneracy towards the bottom of the wells and the 
equidistant spacing achieved high above the barrier top [19|.) 
Now, for the present potential V- W (q) = q 4 — wq 2 , 



y/wq 2 -q 4 dq= -w 3/2 , (79) 



giving \i = | as the growth order for this spectrum. 
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3.2.2 Asymptotic properties 

An essential calculation afforded by the present formalism is the asymptotic eval- 
uation of the spectral determinants for v — > oo. The straightforward behaviour 
lf7U|) of the individual eigenvalues only suggests that D ± (X; v) should closely re- 
late to D^(X\y/v) = det(-d 2 /d<? 2 + vq 2 + A) ± (the operator -d 2 /dg 2 + vq 2 
being equivalent to \fvH2)- We will therefore need the full expressions of 
these harmonic spectral determinants for all v, which follow for instance from 

^■ fm^ 1 <80) 



4 

We will now connect all the determinants through the respective canonical 
recessive solutions <|39[) . which can be fully evaluated (at A = 0) in the WKB 
approximation: "04,a=o f° r the quartic V v (q) = q 4 + vq 2 on the one hand, for 
which n (<7) = (<7 4 + vq 2 ) 1 / 2 and 

/+°° 1 n 3 71 1 

Il ( q )d q =-±(q 2 +v) 3 / 2 ~-?L-l q + 0(p for q^+oo (81) 

^4,o~ (q A + V)- 1/4 e-(9 2 +") 3/2 / 3 for (g 4 + vq 2 ) - +oo, (82) 

and "02,a=o for the harmonic potential (v<? 2 ) on the other hand, for which 
Ilo(g) = w 1//2 (? and 

/ U () (q)dq = -^q 2 2 ,o^K 2 )- 1/4 e"^ 2 / 2 for(V)-+oo. 

(83) 

(In both cases, the large-g behaviour was checked against ea. (|33|l .) 

Now comes the central feature when v — ► oo: while the quartic recessive 
solution "04, o(<7) (initially specified when q — ► +00) has to match a harmonic 
recessive solution upon penetrating the intermediate region 1 -C (7 -C the 
normalization of 7/14,0 (canonically set for g 3> -\A0 feed not match that of "02,0 
(canonically set for 1 <C q and no q* term). A crucial quantity is actually the 
ratio (V^o/V^o), and it simply emerges by reexpanding ea. (|82fl for ^/v ^> q, as 

^4,0 ~ (vq 2 )- 1/4 e-^/s-^/a „ e"^ 2 / 3 V>2,o for 1 « q « ^ (84) 

Thereupon, invoking ea. (|47J) once for the quartic case, then for the harmonic 
case at A = 0, and finally ea. (|80|l . the latest result translates to 

B ± {0;v)~e-* / '/ 3 D£(0\^~e^ / '' 3 ^ 1 '*B£(0), i.e., (85) 

Qi>)~^+ 1/8 e-^ 3 , Qi-(,)~^-/B e -^/3 (7J _ >oo) . 

(86) 
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So, we obtained asymptotic behaviours for these new functions which strongly 
resemble those of the Airy functions (— Ai' and Ai respectively), with an identi- 
cal growth order // = 3/2. Similar reasonings should extend ea. (|86[) to complex 
v with | argw| < tt, and to 



r(i/4) 



w +1 ^ 8 cos 



w 



3/2 



2^ 

r(3/4) 



w 



-1/8, 



3/2 



3 

(87) 
for w — > +oo. 

The asymptotic formula (|78H79f) for the zeros (— Wk) can now be consistently 
regained: either from the asymptotic formula (|87|) on the negative real axis, or 
from the exact quantization condition l|76(l asymptotically expanded by means 
of ea. l|86|l on the half- line {aigv = 7r/3}. 

Fig. 1 plots the pair of functions Qi and their asymptotic forms. Given 
their variegated properties, the idea that Qi 1 *" should be reducible to simpler 
known functions seems unlikely. 



3.2.3 Spectral functions of the zeros and applications 

We must also consider spectral functions of the generalized spectrum {wk}- its 
zeta functions Z ± (s,v), and spectral determinants T> ± (v), which now refer to a 
singular operator, i.e., 

Z ± {s,v) = Tr[(q- 1 H 4 q- 1 +v) ± ]- s , V ± (v) = det {q~ 1 Hiq' 1 + v)* , (88) 

whereas Qi (v) = det(i?4 + vq 2 ^; hence we cannot readily assert that this 
spectrum {wk} is admissible in the sense of Sec.l, nevertheless all the ensuing 
consequences are numerically verifiable and support such an assumption. (If 
this singular operator q~ 1 H4q~ 1 is self-adjoint, then the generalized spectrum 
{u>k} is real as we assumed.) 

Firstly, the semiclassical quantization conditions l|78H79|) for this spectrum 
fix the leading trace identities according to ea. (|36|l : 

6±= T l/4 =>. Z^O) =±1/8. (89) 

Then, even though Qi ± (w) are spectral determinants (at a frozen energy), 
they do not have to coincide with "D^^v). Simply, both being entire functions 
with the same order | and the same zeros, they must be related as T> ± (v) = 
C ± e c±v Qi ± {v) (C±, c± constants). 

Quite generally, a complete identification of such free constants can be based 
on a principle of "semiclassical compliance" whenever the spectrum of zeros is 
admissible: the spectral determinant is a priori known up to a factor exp P(A) 
(with P a polynomial of degree < /i = the growth order) , but only for one such 
P can the canonical semiclassical form JSJl be also satisfied, and this lifts the 
ambiguity completely. 

Here, the asymptotic formulae (|85|l are known for logQi ± (w) down to the 
constant terms (included), and only the latter are actually non-canonical in 
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v — ► +oo, hence necessarily 

V ± (v) = Qi ± (v)/Df(0). (90) 

Several results follow from ea. (|9U[) (and J73J): 

• Explicit Stirling constants for the spectrum {u>fc} (playing the same role as 
V27T for the integers, in view of eq. l|38l) at A = 0): 

V ± (0) ee exp[-Z ± (0)] = Df(0)/Df(Q) (see Table 1). (91) 

• Upon a simple explicit rescaling, the expansions (|13[) written for D (v) yield 

Qi ± ( V ) = Qi ± (0)ex P [-V ^{-v) n ] (\v\<w ). (92) 

n— 1 



Eas. (|91l92|) also amount to 

L-logQi j (0) - I (z±) , (0) + (z ± y(0) (n = Q) (93) 

With these results, we can now draw further consequences from the func- 
tional relation l|74l) and from its coincidence with the homogeneous quartic 
ea. (|143|) up to the exchange of parities. 

• The complete determinants D = D + D~ being unaffected by this interchange, 
the quartic cocycle identity (|145|l of D± must remain satisfied by its analog, 
namely the new product function (Qi + Qi _ ). Hence that functional equation 
11451) now shows two completely different entire solutions, the former having 
order /X4 = | and the latter \i\ = | (tied by the duality relation /xj 1 = 2). 

(We do not know if still other nontrivial entire solutions may exist.) 

• Because of eas. H90fl9"2*)) . it is now the spectral zeta functions Z ± (s) which 
inherit sum rules for s = n > 1, and these have to be the quartic rules ()146|l 
with even/odd arguments swapped, giving 

Z"(1)-2Z + (1) = 

2Z~(2)-Z+(2) = 3[Z-(1)-Z+(1)] 2 (94) 
Z(3) = Z(l) 3 /6-Z(l)Z(2)/2 (etc.) 

(every identity of order 3n expresses Z(3n) in terms of the lower Z(m) in exactly 
the same form as for the quartic zeta value Z^(3n), cf. ea. (|146|l 'l. 

We finally evaluate Z ± (l) in closed form, by analogy with the derivation 
of eg- lim for Z±(l). The values Z±(l) themselves are regularized quantities, 
but the first sum rule l|94|) also implies — Z + (l) = Z p (l), and the latter has 
the (semi)convergent defining series X)fcLo(~l) fe /' u; fc- This series can actually 
be summed when the Wk are more generally defined (for N > 3) as the roots of 
det(— d 2 /dg 2 + \q\ N — wq 2 ) = 0, i.e., as the eigenvalues of the singular operator 
q~ 1 H N q~ 1 (cf. eq.JHBll); hence 

Z£(l) = Tr Piq^HNq- 1 )- 1 = Tr PqH^q (P = parity operator); (95) 
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then, thanks to the explicit formulae I|135I137[) giving the kernel of Hj^ 1 , and in 
full parallel with ea. i|139[l . this turns into the explicit integral 



Z P(1) = ^sin^Tr / [KA2vq 1+N/2 )] 2 q 3 dq, 
Jo 



N + 2 

i.e., a (convergent) Weber-Schafheitlin integral, which finally gives 



(96) 



P sim/Tr 2 _ 8y r(3i/)r(4»/)r(5i/) 
Zn[1 > ~ ~2^T {Zv) r(4^ + i/2) ■ (97j 

For N — 4, by the first of eas. (|94|l . this also yields the special values — Z + (l) = 
-\Z~(l) as given in Table 1 (as well as = -(log Qi+)'(0) = -|(logQr)'(Q), 
thanks to ea. (|9^|) for n = 1). 

Table 1 includes some values of these zeta functions Z ± (s). We computed 
them both analytically and by brute force (using the data from eas. 1(751781 - 
179(1 within eas. ((35138(1 ^. and were thus able to numerically check all the above 
results. 



3.3 A mirror problem without a solution 

The new functions Qi algebraically resemble the spectral determinants of the 
homogeneous quartic potential, but share many qualitative and asymptotic 
properties with the Airy functions, especially their order fix = |. This co- 
incidence clearly reflects the N = 1 N = 4 duality. It is then tempting to 
seek a fourth pair of functions in the empty symmetrical position: verifying the 
Airy functional identities (with +/— arguments swapped), but having the order 
fi4 = | and the qualitative features of the homogeneous quartic determinants. 

However, not only is it difficult to conceive such functions around the spectral 
framework of the linear potential, but in fact it is easily shown that a perfect 
mirror solution cannot exist. The functions of such a pair should have only 
negative zeros (— £fc) to fully resemble the determinants . They would also 
obey the Airy sum rules (|152|l with the +/— superscripts exchanged, i.e., 

Z~(1)=0, Z+(2) = Z+(1) 2 , (etc.), (98) 

with Z ± (s) = J2k £ k S (running over ^f? k) now also for s — 1 since /i = | < 1 
(contrary to the Airy case, where Z ± (l) underwent regularization) . Under this 
precise circumstance, each of the leading putative sum rules above is already 
broken. Therefore, we can get no other partner functions to Qi 1 * 1 in the precise 
manner described. 

(By contrast, this question remains open if we relax the 6k to allow negative 
or complex values, or the order constraint \x = |.) 

3.4 Concluding remarks 

Our study of the functions Qi has remained introductory. We have explicitly 
proved neither that their zeros are purely real (negative), nor that they form an 
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admissible sequence. It would also be nice to know their asymptotic expansions 
I|86I87|) to all orders in v as for the Airy function (those would also give higher 
trace identities for Z ± (— m)); and accessorily, to find closed forms for Z ± (2) like 
those for Z N (2) (JT] App.C). (From a general standpoint, one needs to extend 
the formalism of Sec.l to operators H which can be singular as in ea. (|88|l . so as 
to encompass the Qi* functions; self-adjointness of these operators should also 
be ascertained.) 

The functions Qr* 1 have truly revealed hybrid features. Their analytical and 
algebraic properties are undoubtedly quartic, while their asymptotic properties 
are close to the Airy case. (Unlike the Airy functions, they do not satisfy any 
obvious differential, or linear-difference, equations.) All in all, their structure 
is very simple and strongly reminiscent of the spectral determinants of homo- 
geneous potentials (cf. Appendix), but overall not reducible to the latter; thus, 
Qr* 1 provide seemingly new solutions to the functional identities governing the 
homogeneous quartic determinants. We therefore hope that they might also find 
some roles in the correspondences recently unraveled between those functional 
equations and exactly solvable models of statistical mechanics or conformal field 
theory [Hj. 

4 Exercise III: Supersymmetric 2 binomial po- 
tentials 

Keeping the same techniques as previously initiated for the functions Qf*, we 
now turn to the zero-energy determinants for some other binomial potentials on 
the half-line {q > 0}: namely, 

V(q) =q N + vq M , with N even and M = f - 1 throughout, (99) 

because the corresponding Schrodinger equation (|16fl has special properties: it is 
solvable at A = (in terms of confluent hypergeometric functions |2nH2])j an d 
for selected values of v it also provides the simplest examples of supersymmetric 2 
and of quasi-exactly solvable systems [21]. Then all calculations may strictly 
follow the previous pattern (referring to Sec. 3 for details), yet some of the final 
results will be quite different. 

For the exact formalism of Sec.l, retaining the notations of ea. (|48|l . these 
potentials enjoy an exclusive symmetry, namely, all their conjugate potentials 
are real: 

V [e] (q) = q N + {-lfvq M for all I (100) 

Moreover, the evaluation of ea. (|25[l gives a special non-zero residue formula for 
the first time: (V(q) + A)- s+1 / 2 - q N(-s+\/2) + t ,(_ s + l/2)q~ 1 - Ns for q -> +oo 
implies 

0-i(s) = v(-s + 1/2) (independent of A, N). (101) 
We henceforth focus upon the restricted determinants Dj^(\ = 0; v), where 

D%(X;v) = det(-d 2 /dq 2 + q N + vq^- 1 + \)± . (102) 
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As with Qi^ before, two immediate results are a pair of special values ex- 
plicitly recoverable from eci. (|136l) . 

D±(0; 0) = det(-d 2 /dg 2 + q N ) ± = D% (0), (103) 

and a main functional relation, drawn from ea. (|50fl ) but now taking a special 
form, due to the particular dependence of the degree M upon N and to ea. (|101|l 
(we recall that <p — j^rg): 

e +i ^ 4 D+ (0; -v)D N (0; v) - c^/ 4 D+(0; v)D N (0; -v) = 2ie+W4 . (104) 
4.1 Zeros of D%(0;v) 

The zeros of D^(0; v) are again the values v — —w^ < for which A = is an 
eigenvalue of Hn + v\q\ M (a generalized spectral problem). But now ea. (|104fl 
is very close to the harmonic functional relation l|156|l . especially considering 
the value it of the rotation angle acting on the spectral variable (and also the 
special right-hand-side phase); then just as ea. p56|l . ea. (|104Jl splits into real 
and imaginary parts and reduces to 

J D+(0;«)^(0;- W ) = 2(sin|)" 1 cos (105) 

which is exactly a (shifted) Gamma-function reflection formula. This spectral 
problem is then exactly solvable like the harmonic case: the zeros of the right- 
hand side, which form one doubly infinite arithmetic progression, must simply 
be dispatched according to their signs towards one or the other factor on the 
left-hand side, resulting in the exact eigenvalue formulae 

TV N 
w 2n = — + (N + 2)n, w 2n+1 = — + 2 + (N + 2)n, i.e., (106) 



1 . N-2 



Wk = k+-± — for k = "'i'i>- (107) 

A^ + 2 2 2(iV + 2) 1 - 3 ' 5 '- v ' 

The spectrum {wk} therefore has growth order unity, irrespective of N. For 
N = 2, the generalized spectral problem restores the standard harmonic oscil- 
lator problem det(— d 2 /dg 2 + q 2 + v) = 0. In reverse, the subsequent results will 
prove perfect generalizations to all even degrees N of the classic harmonic os- 
cillator properties. This clearly provides another view on the partial solvability 
properties of the potentials 

For instance, we can show that semiclassical quantization is exact for this 
zero-energy generalized spectrum. We proceed just as for Qi in Sec.3.2.1; for 
large v — —w < 0, now the comparison equation is (— d 2 /dq 2 — wq )if> = 0, 
solved by 

— + 1 i 

Mq)^Q 1/2 J^u(V^^t) (r d d solutions), ^ d -^, (108) 



2 



23 



hence the asymptotic matching in the intermediate region 1 -C q <C w 2v yields 
the semiclassical quantization formula 

while the left-hand side, under the change of variables wx = q~ + , yields 

thus the resulting Bohr-Sommerfeld rule manifestly coincides with the exact 
one GUI). 

The exact quantization formula (|1()7|) — like our earlier eq. I|76|) for the zeros 
of Qi ± (w) = det(— d 2 /dq 2 + q A + vq 2 ^ — displays a hybrid character. The 
right-hand side of ea. (|107|) is that of the exact quantization formula H51|l for 
an arbitrary potential of degree N (whose semiclassical quantization formula is 
definitely different); whereas its left-hand side is linear in the spectral variable 
(entirely coming from the term (— <^/3_i(0)/7r) in eq.JSTJ), and its semiclassical 
form is exact, both as in the harmonic case. 



4.2 Asymptotic properties 

We can get the asymptotic v — ► +oo form of £>^(0;t>) just as in Sec. 3. 2. 2: 
we relate these determinants to the canonical recessive solution $n for the 
potential V(q) = q N + vq M at A = 0, then match '5 at with the solution ipM of 
the comparison potential vq M . Hence we need the canonical WKB forms <|39(1 
for those two potentials, and primarily the symbolic integral V(q') 1 ^ 2 dq' . 

• Full potential V: with a change of variables as above, J V{q) x / 2 dq = 
4i/v J Vl + x 2 dx, we have the obvious primitive 



V{q') l/2 dq' = 2v 



q 



2v 



v arcsinh 
N 1 



,1/2' 



W (^T + 2^° g9 ~ 2 l0gW + 1 ° g2 ) +(?f+1 + 



(q 



we must then set J^°° V(q') 1/2 dq' d = - V{q') 1/2 dq' + C , with the target that 
the large-g behaviour of 'J/jy ~ V(q)^ 1 ^ 4: expf/^ 00 V(q') 1 ^ 2 dq'~\ should obey the 
canonical eq.lgOl), now with e c = 2~ v/N by eas. (1331101(1 : i.e., 



V N {q) 



-v/N -N/i-v/2 -2uq 1+N/2 



(112) 



This adjustment for C according to eas. pillll2|l yields the full WKB specifi- 
cation of the canonical recessive solution, when V(q) = q N + vq M , as 



(111) 

+oo); 



*jv(g) 



-v/N p -i/t)(log'U-l-21og2) 



V(q) 



-1/4 



exp 



Viq'f^dq' 



(V(q) - +00). 
(113) 
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• Comparison with the potential vq M : the corresponding canonical recessive 
solution tpM is expressible from ea. l|135|l exactly, but its WKB form suffices here: 

Vm(?) ~ K M )- 1 / 4 exp[-^f+ 1 /(f + 1)]. (114) 

Now, the matching with eas. (|113llll|) reexpanded for » » 9 yields the result 

V N (q) ~ 2~ v / N e-» v ( lo z v ~ 1 ~ 2lo z 2 1^ M {q) (v -» +00). (115) 

This then translates back to the determinants, by ea. (|47|l . as 

de t(-^ + g Ar + V / ) ± ~ 2-»/*e- w (iog«-i-2io g2 ) det( __ ^, ^m )±> (U6) 
dq z dq z 



But det(-d 2 /d<? 2 + ug M ) ± = v A " z U°) D±(0) by eq-JTSJl, so that finally, 

D±(0;v) ~ e- 1 '^" 1 ' 2^T^% ± ^L>^_ 1 (0), v -► +00. (117) 



4.3 Spectral functions of the zeros and applications 

As in Sec. 3. 2. 3, we also need the spectral determinants T> N {v) built directly for 
the generalized (even and odd) spectra {wk}- Since these form the exact (semi- 
infinite) arithmetic progressions l|l(J7fl . of growth order 1, the answer must now 
have the exact form V^{v) = C ± e c±v /T{v{v T 1) + 1/2)- Compliance with 
eq.® then fixes the constants, giving 



^" 5 rM. T D + i/2) < 118 » 

(which agrees with the harmonic spectral determinants (|155|l for N = 2). 

We can now complete the evaluation of the former determinants D^(0;v) 
themselves, which must have the same zeros as V^(v) and the same order 
unity, hence necessarily D^(0;v) = C± e c±v T>j^(v). This plus the semiclassical 
constraint JHJ again fix the constants, yielding 

D±(0;v) = 2T*WTv v Di (0)V%(v). (119) 

2 

Contrary to ea. (|90|l for the Qi* case, here the latter determinants are known 
by ea. (|118fl . so we end up with fully closed forms, 

2 -v/N^ v y { v+l)+l/2 T (_ 2v) 2 ^/N {Al/r (v-l) + l/2 r{2u) 

^ (0;U) = r(u(v 1) + 1/2) ' ^ (0;U) = !>(« + !) + 1/2) ' 

(120) 

Remark: at v = 0, the identity 1)1191) specifying the ratios C± e c±v — -D^(0; v)/T>^(v) 
simply becomes the duplication formula for T(2v), a fact which also directly fixes 
the constants C±; by contrast, we see no "cheap" way to obtain the constants 
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c± — i.e., c + = c_ = (log2)(2f — = (log2)( 2( ^ + 2 2 ^ ); in particular, ea. (|105|l 
alone is of no avail in this respect. 

Alternatively, we can invoke the solvability of this potential at zero energy 
to obtain the spectral determinants directly, just as we deduced ea. ((136|) for the 
homogeneous potentials (cf. Appendix). Here the canonical recessive solution 
\& jv normalized by eq. 1(112(1 can be expressed indifferently in terms of a confluent 
hypergeometric function U(a, b, z) or a Whittaker function W K ^(z), as ([23] eq. 
2.273(12)) 

* N (q) = 2- v ' N {AuY^ +l ' 2 e- 2 »i 1+N/2 U{u(v - 1) + 1/2, 1 - 2u, W +N ' 2 ) 
= 2- v / N (4vrq- N / i W„ uv ,„(W +N/2 ) (121) 

(the normalization is fixed by reference to the known q — > +oo forms of these 
functions). In turn, the connection formula U(a, b, z) = T(l — b)M(a, b, z)/ r(l + 
a - b) + z 1 ~ b T(b- 1)M(1 + a - b, 2 - b, z)/T(a), together with M(a, 6,0) = 1, 
determine the values of \1/Ar(0), '^(O), which finally yield the determinants 
through ea. 1)47(1 : these then coincide with ea. 1)120(1 indeed. (Conversely, the 
present method can be viewed as a purely spectral derivation of the connection 
formula for that confluent hypergeometric function. Again, however, the right 
reference normalization 11121) of <3>Ar has to be expressly fed in.) 

(Remark: the knowledge of ^N(q) also yields an expression for XX - ^) k / w k, 
in full analogy with eas. l(95hl^7|l above, but here this only recovers a special 
case of the known integrals J °° W KiV (z) 2 dz/ z ([2], eq.7. 611(4)); likewise, the 
sum rules for ^2(±l) k /w% (n > 1) should only reproduce elementary identities 
among Hurwitz zeta-values here.) 

For N = 2, eq. l(12(J(l restores the ordinary harmonic spectral determinants 
11551) . The next even case, V(q) = q 6 + vq 2 , is strongly highlighted in the studies 
on quasi-exactly solvable potentials |22j . However, our present results hold 
identically irrespective of the parity of the full potential q N + vq M , beginning 
with the potential q 4 + vq on the half-line {q > 0}, so we end with a few remarks 
about these non-even potentials. 

4.4 The case of non-even potentials 

Our exact results above hold equally well for non-even potentials V(q) — q N + 

N -i 

vq~~ , obtained when TV is a multiple of 4, as for even ones. Then, as always, 
the exact zeros {wk} of det(— d 2 /dg 2 + q N +vq~~ 1 ) ± as given by ea. 1(107(1 refer 
to the potential defined on the half-line {q > 0}, with a Neumann/Dirichlct 
condition at q = for the +/— parity, or equivalently to the singular potential 
y(|g|) over the whole real line. 

For N multiple of 4, an additional exact spectral property is derivable. We 
now consider the (non-even) potential V(q) over the whole real line. The com- 
plete spectral determinant D for such a potential of even degree N is not given 
by D(X) = D + (\)D~(\) (which corresponds to the even potential but 
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in full generality by ^U] 



D(A) = ^[D+(\)D^ 1+N ^-(X) + D^ 1+N ^+(X)D-(X)}. (122) 

In particular, for V(q) — q N + vq^^ 1 with N = [mod 4], by ea. 1(100(1 . 

I cos TTiyv 

B N (0:, v) = - [D+(0; v)D^(0; -v) + D+ (0; -v)D N (0; v)} = - (123) 

z sin ttv 

(the explicit end result uses ea. l|12()[l ). Thus, the values of w = — v such that the 
potential q N — wq^^ 1 on the whole real line has a zero mode are also exactly 
quantized: they are the zeros of cos N \ 2 V > ^ e -> 

w' n = {N + 2)(n + 1/2) (ti GZ) if N = [mod 4], (124) 

(a spectrum naturally invariant under reflection; whereas in the even potential 
case N = 2 [mod 4], the spectrum {wk} of eq. l|107|l is recovered). 

Curiously, the eigenfunctions 1(121(1 corresponding to the explicit generalized 
eigenvalues 1(124(1 do not seem to reduce to elementary functions, while they do 
so for the generalized eigenvalues ((107(1 : still, the basis that they form might 
prove useful (e.g., the N = 4 basis for general quartic anharmonic oscillators). 

In conclusion, a unified analytical formalism has displayed the quantization 
of the special zero energy in potentials of the form q N + vq~~ x (beginning 
with TV = 4), as well as the connection formula for the corresponding confluent 
hypergeometric functions, to be clear generalizations of the explicit exact quan- 
tization scheme for the harmonic oscillator. Contrary to the previous example, 
however, this analysis treats the problem entirely by known functions and does 
not generate any new ones. 

We are grateful to A. Turbiner for suggesting us to apply exact WKB analysis 
to quasi-exactly solvable potentials. 



A Appendix: formulae for the potentials q 

We recapitulate specific results and formulae, otherwise scattered in many ref- 
erences El Ell El 13 EH IIH El E3 E3 EH1 , about the spectral functions of the 
homogeneous Schrodinger operators on the half-line, 

~ h f d 2 

H N = -—r + q N , qe [0, +oo) N > 1 integer. (125) 
aq z 

(Not all iV-dependences will be systematically stated.) As in Sec. 1.2, we call 
H^f (resp. H^) the operator with the Neumann (resp. Dirichlet) condition at 
q = 0. 

Three cases of special interest to us will be, on one side, the harmonic (N = 
2) and linear (N — 1) cases, both of which are describable using known special 
functions, as opposed to the non-elementary quartic case (N = 4). At the same 
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time, N — 4 and N = 1 are both regular cases, and dual to each other (they 
share the same number of conjugates L = 3), whereas N = 2 stands out as a 
singular (confluent) case (also self-dual) 



A.l General TV 

Notations: 



J D±(A) d ^ f det( J ff±+A), and ^ = ^±H, y, = , (I2(,i 



The set of exponents in eq.@ reduces to {0} U {(2n— l)/i} n= o,i.2,- -; the coeffi- 
cient of the leading singularity t _Al is 

c% = (2V^)- 1 T(1 + 1/N), (127) 

In accordance with eq. (J5J , the trace identities then have the pattern [5J |3] 

^(0) = 0, Zjr(0) = 1/2 (128) 



Z N (-m) = unless m = (A + r)(l + |r 
i?^(— m) = unless m = r(l + 



for m and r G N (129) 



• As all the conjugate potentials V^(q) (eq.l^HJ) coincide here, the main 
functional relation l|5U|) for the spectral determinants boils down to 2 

e+ i ^ 4 £»+(e-^A)£»^(A)-e- i ^ 4 J D+(A) J D^(e- i ^A) = 2i (JV ^ 2) (130) 

or equivalently, to a multiplicative "coboundary identity" linking the full and 
skew determinants, [23 El EH] 

£>£(A)/£>£(e-^A) =e i (-2*^(A)+W2) 5 $Jv(A) = f arcsin [D N (e~ iv X)D N (X)]~ 1 / 2 , 

(131) 

the branch of the arcsin being fixed at A = with the help of eq. I|13UI) : 

D N (0) = (sin ip/4)- 1 =► * JV (0) d = (132) 



The exact quantization condition, drawn from ea. l|13U|l . is |271 128ll2*^| 



2 1 N — 2 

-arg£>±(-e- i ^) B=Efe =fc + -±— _ for fc = (iV ? 2). 

(133) 

• By multiplying together all the conjugates of ea. l|131|l . a corresponding 
multiplicative "cocycle identity" (of length L given by ea. (|49|l ) results, providing 
a consistency condition upon D alone in the implicit form |26| 

L-l 

^^(e-^A)=^/4, (134) 
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i.e. an autonomous functional equation, circularly symmetric of order L (and 
convertible to a polynomial form) for the complete determinant Djy. 

• Some special values of spectral functions also become explicit d] 
App.C) thanks to the solvability of the differential equation (H^ + \)ip\(q) = 
at A = 0, specifically here by Bessel functions (^\, vol.11 ch.7.2.8): the canonical 
recessive solution l|3l3|) . of asymptotic behaviour ijjo(q) ~ q~ N / i exp[— q 1+N / 2 / (1+ 
AT/2)] prescribed by eq.QU, is 

MQ) = 2^d 1/2 KA.W +N/2 ) {v = j^) (135) 

(normalized by reference to the known z — > +oo behaviour of K u (z)). In turn, 
the values V'o(O), V'd(O) follow from K v {z) — 7r(2sin vk)^ 1 [/_„ (z) — I v (z)] plus 
I± v {z) ~ {z/2) ±u /T{±v + 1) for z -> 0, and finally eq.gTJ) yields 

^(0)= ^Ll^- ^(0) = F(iy) ^ /2 ^(0) = ^— (136) 

V Nv l Z ^J-K \/TT Sin 

(the last result agreeing with ea. l|132|) '). 

Ea. l|135|l also generates formulae for Z^(n), n = 1, 2, • • •, by fully specifying 
a general expression of the ID Green's function, 

(qKH + Xy'W) = W r A - 1 VA(-min{ (? , (Z '})^ A (max{g,q'}) (137) 

(W x = Wronskian{^ A (- (7 ),VA (<?)}) (138) 

when H = H N and A = (with W\ = 2D(X) by eq.lgSJl, implying W = 
2(sin <p/4) _1 by ea. (|136|> ). Integral formulae giving Zm (n) as Tr^ 1 )", and 

Zjf(ri) as Tr P{P[J j 1 ) n (P d = the parity operator), thus become explicit. For 
n = l, the latter is the simpler: 

A /"CO -| 

Z P(1) = -sin^TT / [K v {2vq l+N / 2 )] 2 q dq, u = ——; (139) 

ft Jo i\ + z 

finally this evaluates by a Weber-Schafhcitlin formula , vol.11 ch.7.14, eq.(36), 
and likewise for Zn(1) 3 ): 

zin) = f^(2.) 2 - r( ; ) , 2 r(2 ; ) 1 ^ ) , ^(D = ^^(i). (ho 

2^^ r(2f + l/2) tan/y7r 

(Remarks: - the results (|136f) . (I140|l also work for N = 1, 2; - this approach 
can still handle Z^(2), but with final results reducible only to 4F3 generalized 
hypergeometric series: |TT], App. C.) 

• For general integer n, by contrast, the farthest explicit algebraic result 
we can reach is a single identity (a sum ruJe) at the level of each doublet 
(Ztf(n), Zjf(n)), and this comes simply by expanding the functional identity 
ifl^njl in all powers A", as [TT] 



sm 



f + £ sin J ^M(-A)"] = exp [ Z ' N (0) + £«r-^ 

n— 1 m— 1 



2 m v ' 

(141) 
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then equating both sides of this generating identity at each order n: at n = 
we already obtained Dn(0) = l/sin(y>/4) ('ea. (|132ll or (SEJ)j then, at higher 
orders, we get sum rules in the form: 

cot | S i n -cos | Z N {1) = (W^2) 

cot | sin (2) -cos ^^(2) = [2 cos | Z P N (l)] 2 (142) 

cot|sin^Z^(3)-cos^Z A r(3) = 4cos |[3 cos | Z&(1)Z&(2) - 2 cos 2 | Z^(l) 3 ] 



cot — sin ^ Z^(n) - cos ^ Zjv(n) = a polynomial of {Z^(m)}i< m< „ 

(the first line also follows from eq.l^UJ)). 

Other results are resurgence properties of the spectral determinants (namely, 
exact analytical constraints on their (l/?i)-Borel transforms) |llll2"ll3"U|. 



A. 2 Special cases 

Out of the preceding results valid for general N, we only restate those which 
are specifically needed in the main text, or which take a particular form; special 

and numerical values are in Table 1. We will be using j = f e 2l7r / 3 . 



A.2.1 N = A 

This case (quartic oscillator) has order (j, = ~. 
Main functional relation (|130|l : 

c +W6^+ (j -i A) ^- (A) _ e -W6^+ (A) ^- (j -i A) = 2i ^ = 2tt/3). (143) 
Exact quantization condition (|133|l : 

~ arg^±(-j 2 £;) B=Bfc = k + \ ± i for k = (144) 

Cocycle functional equation l|134|) : its algebraic form is 

£ 4 (A).D 4 (jA).D 4 (j 2 A) = D 4 (A) + D 4 (jA) + Z? 4 (j 2 A) + 2. (145) 

Special values for Df(0), Zf(l): see Table 1 (note that D 4 (0) = 2). 
The sum rules l|142H can also be written for N = 4 as 

Z+(1)-2Z 4 -(1) = 

2Z+(2)-Z 4 "(2) = 3[Z+(1)-Z 4 -(1)] 2 (146) 
Z 4 (3) = Z 4 (l) 3 /6-Z 4 (l)Z 4 (2)/2 (etc.). 

(In addition, Z^(s) have been asymptotically evaluated for s — ► — oo, by 
exploiting the resurgence equations for the spectral determinants. 012) 
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The eigenvalues E k are the (unsigned) zeros of the Airy functions, here meant 
as Ai (for odd parity) and Ai' (for even parity) ( 21 , ch.10.4). A few results 
nevertheless seem new [221 • The order is /j, = |. 

The canonical recessive solution simply relates to the Airy function: 

^ x {q)=2^M{q + X). (147) 

The spectral determinants are then also Airy functions by ea. (|47|l . 

D+{\) = -2V^Ai'(A), L»r (A) = 2V^Ai(A) => D X (X) = -27r(Ai 2 )'(A), 

(148) 

and their functional relation (|130[) boils down to the classic Wronskian identity 
for Ai(-) and Ai(j 2 -): 

c+ i -/ 3 ^+(jA) J Dr(A)-c- i ^ 3 ^+(A)^r(jA)^2i {<p = 4rr/3). (149) 

Exact quantization condition (|133f) : 

\r g D±{-]E) E=Ek = k + ^ T \ for * = Sfc: (150) 

Cocycle identity: the algebraic form of (|134|l for N = 1 is 

^ 1 (A) 2 +I? 1 (jA) 2 +C 1 (j 2 A) 2 -2[ J D 1 (jA) J D 1 (j 2 A)+^ 1 (j 2 A) J D 1 (A)+ J D 1 (A) J D 1 (jA)]+4 = 

(151) 

Eas. (|136| ) for Df (0) just express the known values of Ai(0), Ai'(0) (see Table 

!)■ 

Sum rules JTHJ: 

Z+(1)=0, Z-(2) = Z-{\)\ Z 1 (3) = 5Z 1 (l) 3 /2-3Z 1 (l)Z 1 (2)/2, (etc.). 

(152) 

Here, due to a supplementary set of relations, all values Z^(n) become recur- 
sively expressible as rational functions of 

r d = Dl (0) = - Ai'(Q)/ Ai(0) = 3 1/3 r(2/3)/ r(l/3) « 0.729011133 ; (153) 

they are all irrational except Z^(\) = (a regularized value, however) and 
Zj~(3) = 1 (Table 1 lists all Z^(n) analytically up to n — 3). 

A.2.3 A = 2 

The finite- N case where the spectrum is known (Ek = 2k + 1), of order fi = 1 
(an integer, which makes it a singular case). 

The canonical recessive solution relates to a parabolic cylinder function ( |2J , 
ch.19; .20 , vol.11 ch.8) 

M<l) = 2^' A U(X/2, V2q). (154) 
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The determinants are given by Gamma functions, 

Da(A) - rm ' 2( )_ T(W ^ 2(A) - r(^) ' 

(155) 

their functional relation is now drawn directly from ea. l|50(l instead of H130[l : 

C+W4^+ ( _ A)jD - (A) _ e -W4^+ (A)C - ( _ A) ^ 2ie i.A/4 ^ = ?r)) (156) 

and its right-hand side carries a special extra factor because for homogeneous 
potentials, 

/3_i(s) = 0, except: 0-i(s) = A(— s + 1/2) for N — 2. (157) 

For this special value ip = tt, the functional relation further splits into its real 
and imaginary parts at real A, and thereby reduces to 

£>+(A)£>a (-A) = 2 cos J(A - 1); (158) 

hence it just amounts to the reflection formula for T(z) (like the cocycle identity, 
which we do not write). 

Here, the exact quantization condition Ek = 2k + 1 trivially comes by dis- 
patching the obvious zeros of the right-hand side of eq. i|158H : the positive ones 
to D^i—X), the negative ones to -D^~(A). 

The anomalous value 1|157|) accounts for special N = 2 trace identities: 

Z 2 ± (0,A) = (-A±l)/4; (159) 

the anomaly vanishes at (and only at) A = 0, and indeed the general formulae 
(14711361) for D%(\ = 0) correctly agree with eq. I|155(l to give 



D 2 (°) = -7TT ~ 0.977741067, D^{0) = -jj- w 1.446409085 =*> D 2 (0) = 

(160) 

The spectral zeta functions are 

oo 

Z 2 (s) = (l-2- s )C(s), Zl{s)^P{s) ( d ^^(-l)V(2fc + l) s ); (161) 



k=0 



their sum rules only reproduce the known special values C(2n) and (3{2n + 1) 
(EU, ch.23) 

Z$(l)=n/4, Z 2 (2)=tt 2 /8, ^(3) = tt 3 /32, (etc.). (162) 
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Figure caption 

The functions Qi~ and (— Qi + ) (of ea. (|72|l L which specially resemble the 
Airy functions Ai and Ai' respectively, plus logarithmic plots to magnify their 
large- u asymptotics on the decaying side (v » 0). Numerical values: + for Qi + , 
and o for Qi~ (computed using ea. i|38|) . with an input of K 10 3 numerical 
eigenvalues of the operator (— d 2 /dq 2 + q 4 + vq 2 ) for every value of the coupling 
constant v). The dashed lines plot the large- v asymptotic formulae I|86I87|) . 
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